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Heavy quark production probes QCD at the interface of the perturbative and non-perturbativc 
regimes. Studying the production of heavy quarks is an important test of models in both 
regimes. In this article, recent results on beauty and charm production from the CDF and 
D0 experiments at the Tevatron are reported. These include measurements of correlated bb 
production, the t/>(2S) production cross section, and T(1S) and T(2S) polarization. 



1 Introduction 



The study of heavy quark production is an important test of models of both perturbative and 
non-perturbative quantum chromodynamics (QCD). Since the development of nonrelativistic 
QCD (NRQCD), the agreement between theoretical predictions and measured values of heavy 
quark production cross sections has improved greatly. However, recent J ftp and ip{2S) polar- 
ization measurements at the Tevatron^ indicate the understanding of heavy quark production 
is not yet complete. 

The Tevatron at the Fermi National Accelerator Laboratory produces pp collisions with a 
center of mass energy of 1.96 TeV. The Collider Detector at Fermilab (CDF) and D0 experiments 
employ general multipurpose detectors to reconstruct pa rticl e physics events from these 
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collisions. With 6-hadron cross sections of ~30 /ib [\rj\ < l.O),^ 1 -^ CDF and D0 have a wealth 
of experimental data on 6-hadrons. 



2 Correlated bb Production 



The cross section for producing, in hadronic collisions, both the b and b quarks centrally and 
above a given pr threshold is referred to as the bb correlation a bb . The exact next-to- leading 
order (NLO) prediction of a bb appears to be a robust perturbative QCD prediction, differing by 
only a few percent from the leading order (LO) prediction. However, Run I measurements of 
a bb at the Tevatron are inconclusive, with an average ratio of the measured a bb over the exact 
NLO prediction (R 2b ) of 1.8 with a 0.8 RMS deviation^ 

In this measurement,^ a bb is obtained in 740 pb -1 of CDF data using dimuon events with 
Pt(m) > 3 GeV/c, \r}\ < 0.7, and invariant mass € [5,80] GeV/c 2 . This corresponds to bb 
pairs with pr > 2 GeV/c and rapidity \y\ < 1.3. At the Tevatron, dimuon events mainly result 
from the decay of heavy quark pairs (bb or cc), the Drell-Yan process, the decay of charmonium 
and bottomonium, and the decay or misidentification of ir or K mesons. To determine the bb 
and cc content of the data, we fit the impact parameter distribution of the muon tracks. The 
impact parameter do is defined as the distance of closest approach of the track to the primary 
event vertex in the transverse plane, and it is proportional to the decay time of the parent 
particle. 

The one-dimensional impact parameter distributions of muons from b and c decays are mod- 
eled by a tuned Herwig simulation,^ while the distributions of prompt muons are reconstructed 
using muons from T(IS') decays in the data. Because the impact parameters of the two muons 
are to first order uncorrelated, the three ID templates may be combined into six 2D templates 
to represent each possible dimuon source (bb, cc, cb, prompt-prompt, prompt-6, and prompt-c). 
These six templates are then used to perform a maximum likelihood fit to the 2D distribution 
of the impact parameter of both muons to extract the bb and cc components. The projection of 
the 2D impact parameter distribution is compared to the fit result in Figure HJ 

From this fit, we measure the dimuon cross sections to be <7 6 = 1549 ± 133 pb and 

<7 c ^ Mi c^ = 624 ± 104 pb, where the quoted error is the sum in quadrature of statistical and 
systematic uncertainties. In order to compare with theoretical predictions, we evaluate the NLO 
dimuon cross section using the MNR generator with the events decayed by EvtGen.^ This gives 
a value of <r^ = 1293 pb, resulting in a ratio of i?26 = 1-2 ± 0.2. This measurement is 
in agreement with the NLO theoretical prediction, and does not confirm the anomalously high 
dimuon cross section observed in Run I. 

3 ip(2S) Production Cross Section 

Charmonium production provides another arena in which to test our understanding of QCD. 
The development of NRQCD was prompted in part by the CDF Run I measurements of J/tp 
and ip(2S) production cross sections.^ A new CDF measurement of the ip(2S) production cross 
section pushes the p^ range farther into the perturbative QCD regime than was possible with 
Run I data.El 

We reconstruct ip(2S) — > fi~*~ fi~ using 1.1 fb _1 of CDF data. We then perform an unbinned 

maximum likelihood fit in the ip(2S) mass and proper decay length ct distributions. The mass 

fit separates the signal from the background, while the ct fit separates promptly-produced ijj(2S) 

from ip(2S) originating from secondary decays of long-lived particles (predominantly B mesons). 

The i/j(2S) pt range of 2 to 30 GeV/c is divided into 25 bins, and the signal yield and prompt 

fraction in each p^ bin are extracted by the likelihood fit. 

The ip(2S) acceptance depends upon the ip(2S) polarization. The CDF measurement of 

m 

tp(2S) polarization is statistically limited, with only three measured data points. 1 ^ We take the 
average of these three data points as the effective polarization a c s = 0.01 ± 0.13, where a is 
defined according to Equation ([1]), and use this value to calculate the ~4>(2S) acceptance in each 



bin of ip(2S) pr- 

The final result is shown in Figure EJ where the promptly produced ip(2S) are separated 
from those produced in B decays. The CDF Run I measurement is also shown in Figure [2j The 
integrated cross section has increased by 18 ± 19% from the Run I measurement, compared to a 
theoretical prediction of 14 ± 8% for the change in center of mass energy from 1.80 TeV to 1.96 
TeV.^1 The integrated inclusive differential cross section is measured to be 

a(pp - ^(2S)X, \y(i>(2S))\ < 0.6, p T > 2GeV/c) %/i=1 . 96TeV ■ Br^{2S) -> /X+/0 

= 3.17 ± 0.04(stat.) ± 0.28(syst.) nb. 
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Figure 1: Projection of the two-dimensional impact pa- 
rameter distribution of muon pairs compared to the result 
of a fit for the bb and cc components. The different letter 
combinations in the legend indicate the dimuon source, 
as explained in the text: "B" for a muon from a b quark, 
"C" for a muon from a c quark, and "P" for a prompt 
muon. The legend gives the number of events found by 
the fit for each dimuon source, along with the statistical 
uncertainty. 
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Figure 2: The differential cross section for prompt tp(2S) 
(solid circles) and ip(2S) from B meson decays (hollow 
circles). The cross section for %j){2S) from B meson de- 
cays is scaled down by a factor of 10 so that the two 
curves may be seen separately. Triangular markers indi- 
cate the results of a CDF Run I measurement, which has 
not been scaled to account for thejlifference in center of 
mass energy.^ 



4 Polarization of T(15) and T(25) 

Recent measurements of the J/^J and ip(2S) polarization at CDF show significant longitudinal 
polarization with increasing pp; while NRQCD predicts transverse polarization at sufficiently 
high pt for S-wave quarkonia produced in pp collisions. Now D0 presents new measurements 
of the T(15) and T(25) polarization, another important test of the theoretical approaches to 
QCD. 

The polarization is measured in the parameter a, defined as: 

a T -2a L 

a = (1) 

where cjt and ox are the transverse and longitudinally polarized components of the cross section 
respectively. If the transverse and longitudinal states are equally populated, one measures a = 0; 
for longitudinal polarization, a < 0, while for transverse polarizaton, a > 0. In quarkonia decay 
to a lepton and anti-lepton, a may be obtained from the angular distribution 



Ul\ 9/1* 

— — -r oc 1 + a COS 

a(cos B*) 



(2) 



where 9* is the angle between the T in the lab frame and the positive lepton in the T rest frame. 

Using 1.3 ftr 1 of D0 data, we find 420 000 T(nS) -> / u + >~ candidates. 9 The T data are 
divided into several bins in T (Pr) an d I cos $*|> an d the number of T(1S) and T(2S) in each 
bin are extracted from a fit to the T mass distribution. The mass signal consists of three peaks, 
the T(15), T(2S), and T(3S), where the mass differences between the peaks are fixed to the 
measured values. Unfortunately, the number of T(3S) was insufficient to extract the angular 
distributions. The angular distribution in each bin is compared to T(1S) and T(2S) Monte 
Carlo samples which were generated with the parameter a set to zero. The measured value of 
a for the data is then determined by reweighting the angular distributions in the Monte Carlo. 
The dependence of a on p^, is plotted in Figure [3]for both T(IS') and T(25), along with various 
theoretical predictions. While statistics for the T(25) are insufficient to draw a conclusion, in 
the T(IS') there is significant dependent longitudinal polarization which is only marginally 
consistent with any of the theoretical predictions. 
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Figure 3: On the left is shown the dependence of the polarization paramete r a on for the T(IS'). The black 
points are the D0 data , w hile the green triangles are the CDF Run I result. The yellow band corresponds to 

ii~2i n~si 

an NRQCD prediction, 1 — 1 while the two dashed curves are two limiting cases from &t factorization models. 1 — ^ 
The lower line corresponds to the quark-spin conservation hypothesis, while the upper line corresponds to the full 
quark-spin depolarization hypothesis. On the right is shown the dependence of a on f° r the T(2S I ). The blue 
points are the D0 data and the yellow band is an NRQCD prediction. 



5 Conclusions 



Recent measurements of the bb correlation and the ip(2S) cross section at the Tevatron are 
in agreement with NLO and NRQCD predictions. However, measurements of quarkonia po- 
larization in the same perturbative regime show discrepancies from theoretical predictions. 
Theoretical models are now challenged to match the polarization measurements while continuing 
to describe the cross section data. 
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